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An  accurate  determination  of  the  electromagnetic  scattering  cross  section  of  many  realistic 
targets  must  account  for  their  being  inulti-diinensional  and  penetrable  by  the  incident  EM 
wave.  An  existing  volume  integral  equation  formulation  of  EM  fields  within  continuous 
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source  region  of  this  formulation  to  be  a  dielectric  body  in  whirh  currents  are  induced  by 
the  Incident  wave,  and  which  then  "re-radlates"  the  scattered  wave.  The  relation 
Eioi  “  E-in,.  +  Efi, yields  a  matrix  equation  to  be  solved  for  the  unknown  currents,  which  in 
turn  directly  determine  the  scattering  cross  section.  Numerical  n  .ailts  are  ])ic.seiitcd  for  a 
homogenous,  two-dimensional,  rectangular  target  for  TE  and  TM  polarizations.  The  role  of 
the  source  dyadic  term  In  improving  the  accuracy  of  the  scattering  cross  section  is  investi¬ 
gated  both  theoretically  and  numerically. 
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9.  Upper  Left:  RCS  of  an  Infinite  Square  Cylinder.  0.05  x  2.50  X.  e/co  =  4.0, 

Z  Polarization.  With  Vertical  Target.  Upper  right;  with  horizontal  target:  for  com¬ 
parison  with  upper  left.  Lower  left:  calculated  scattering  patterns  of  a  homogeneous 
plane  dielectric  slab  of  the  same  thickness  and  width  as  above,  with  a  plane  wave 
having  normal  incidence,  and  lower  right:  with  a  plane  wave  at  grazing  incidence 

10.  Same  as  Figure  9,  but  for  xy,  or  TE  Polarization 
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Electromagnetic  Scattering  From  Dielectrics  - 
a  Two-Dimensional  Integral  Equation  Solution 


1.  INTRODUCTION 


There  has  been  considerable  interest  througljoiit  the  1970s  and  1980s  in  the  scattering  of 
electromagnetic  waves  from  dielectrics,  for  a  wide  variety  of  applications.  These  include 
propagation  through  rain  and  snow,  ilctection  of  airborne  particulates,  coupling  of  missiles 
with  plasma  plumes  or  dielectric-filled  apertures,  performance  of  communications  antennas 
in  the  presence  of  dielectric  and  magnetic  inhomogeneities,  scattering  from  birds  and  other 
interference  in  radar  systems,  microwave  hazards  to  biological  (e.sjx’cially  humanj  tissue, 
hypothermia  treatment  of  cancerous  tissues,  and  last  but  not  least,  detection  ol'  non-metaltic 
(that  is.  electrically  penetrable)  components  of  military  aircraft,  missiles,  and  radar  antennas. 
Numerous  papers  have  addressed  mostly  differential,  but  also  some  integral  ctpiation 
techniques  for  calculating  electromagnetic  fields  scattered  by  inhomogeneous  dielectric 
objects.  The  volume  integral  equation  method  is  a  viable  alternative  to  differential  methods, 
but  it  has  not  usually  been  ai)plied  to  multi-wavelength  objects  because  of  the  laigi-  matrix 
inversions  required.  However,  with  the  recent  availability  of  parallel  comimters  and  new 
matrix-solving  algorithms,  this  limitation  is  being  rapidly  overcome 

A  further  limitation  of  integral  equations  concerns  the  careful  treatment  of  the  source 
region,  a  controversial  issue  in  the  recent  past.  Originally,  this  report  was  intended  as  a 
continuation  of  the  source  region  problem  -  atr  application  meant  to  ilemonstrate  a  difference 
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when  one  correctly  represents  the  coincidence  of  the  source  and  observation  points  in  forming  an 
expression  for  the  scattered  electric  field.  The  derivation  of  the  matrix  equation  presented  here 
and  the  computer  program  used  for  the  calculations  thus  follow  directly  from  previous  work  on 
the  source  region  problem.' 

Upon  reviewing  the  existing  literature  about  scattering  from  dielectrics,  one  quickly 
discovers  that  in  almost  every  article  the  primary  reference  is  to  Richmond's  work  of  the  early 
1 960s. *  Since  Richmond's  results  match  the  physical  optics  approximation  results,  they  have 
been  used  as  a  baseline  for  comparison  in  all  subsequent  work.  Surprisingly  enough,  the 
theoretical  development  in  this  report  parallels  Richmond's,  produced  long  before  the  source 
region  controversy  emerged.  Also,  the  results  of  the  computer  program  precisely  match 
Richmond's  results  for  the  two-dimensional  problem.  Since  Richmond  omits  nearly  all  of  the 
details  of  his  calculation,  it  is  difficult  to  determine  whether  he  actually  did  consider  the  source 
region,  or  simply  had  keen  enough  physical  intuition  to  choose  correct  forms  of  the  expressions 
based  on  knowledge  of  their  behavior  for  limiting  cases.  In  this  report,  the  details  of  the  correct 
closed  form  integral  equation  solution  are  presented  explicitly,  and  the  code  is  available  (to  DoD 
workers)  as  an  in-house  resource  for  future  modification,  if  desired. 


2.  OUTLINE  OF  THEORY 


Figure  1  shows  the  geometiy  of  a  bistatic  radar  cross  section  problem.  An  incident  wave 
excites  polarization  currents  in  the  dielectric  target  which  then  re-radiates  the  energy  in  the  form 
of  a  scattered'  field.  Thus,  the  target  is  considered  the  source  of  the  scattered  field,  and  Llie 
standard  cqiiation  Es<.,t  =  E,,„  -  Ej,,,.  is  used  to  calculate  this  scattered  electric  field  at  any  point  in 
space.  When  E^,  ;,!  is  then  expressed  in  terms  of  the  unknown  induced  polarization  current,  it 
(linckly  yields  the  scattering  cross  section,  or  RCS.  used  to  characterize  the  target  s  visibility  to 
radar. 

An  exact  expression  for  the  scattered  electric  field  due  to  currents  induced  in  a  dielectric  target 
has  been  previously  derived,  taking  care  to  properly  account  for  the  singularity  occurring  in  the 
electric  Green's  function  when  the  source  and  observation  points  coincide  (at 
r  =  r'). 


'  Silherstein.  Mariiui  (1991)  Application  of  a  geticralized  Leibniz  ride  for  calculating 
electromagnetic  fields  within  continuous  source  regions.  Radio  Science.  26  (No.  1 ):  183- 190 
(a  more  detailed  version  is  contained  in  Rome  Air  Development  Center  Technical  Report  88-333. 
ADA2  12470. 

^  Richmond.  J.H.  (1975)  Scattering  by  a  dielectric  cylinder  of  arbitrary  cross-section  shape, 
IEEE  Transactions  on  Antennas  and  Propagation.  AP-13;334. 

*  Richmond.  J.H.  (1966)  TE-wave  scattering  by  a  dielectric  cylinder  of  arbitrary  cross-section 
shape.  IEEE  Transactions  on  Antennas  and  Propagation.  AP-14:460. 


2 


+  X 


of  a  Blstatic  Scattering  Experiment 


Ej,j,.,,(r)  =  i(i)Ho  lini  f  G,^  (r,r')  ■  J(r')  dV 
6—0  Jv-Vj 


Lj  •  J(r) 

+  - 

ilOEo 


(1) 


This  outcome  of  the  controversial  "source  region  problem"  contains  the  important  source 
dyadic,  or  L^.  term,  accounting  for  the  field  contribution  at  r  =  r'.  The  calculation  in  this 
report  starts  from  this  expression,  and  solves  the  integral  equation  resulting  upon  its 
insertion  into  the  simple  scattered  field  equation  above.  The  reader  is  urged  to  review 
Yaghjian's  or  Silberstein's  development  of  Eq.  {1).*  '* 

The  following  derivation  is  performed  for  a  two-dimensional  target,  that  is,  an  infinite 
cylinder  of  arbitrary  geometrical  cross  section.  Both  TE  and  TM  polarizations  are  addressed, 
with  plane  wave  incidence  assumed.  The  three-dimensional  problem  would  proceed  in  the 
same  manner,  but  with  a  simpler  Green's  function,  and  is  discussed  briefly. 


3.  THEORY  FOR  TM  POLARIZATION 


We  begin  by  considering  the  details  of  Eq.  (1).  which  can  be  transformed  via  the  Helmholtz 
wave  equation  into  the  form: 


E  (r)  =  1 -  lim  f  [V'xV'x  J(r')]ij)(r.r’)  dV 

6-0  Jv-v, 


L6  -  J 

iuuo 


(2) 


Note  that  the  source  dyadic  term  accounts  for  the  contribution  to  the  scattered  electric  field 
from  the  principal  volume.  V6,  excluding  from  the  integral  the  singidarity  caused  at  the  source 

point.  Also  note  the  dyadic  L,^  varies  with  the  shape  of  this  volume  and  the  fixed  liosition  of 
the  singularity  within  it.  Because  we  have  a  thr<-e-dimen.sional  volume  integi.il  c(|uatir)n  but 
wish  to  investigate  a  two-dimensional  target  as  the  source  of  the  scattererl  field,  we  first 
eliminate  the  z-dependence  by  performing  the  z-integration  of  the  Green  s  function,  and 
expressing  the  current  and  its  derivative  operators  as  two-dimensional  vectors.^ 


■*  Yaghjian,  A.D.  (1985)  Maxwellian  and  cavity  electromagnetic  fields  within  continuous 
sources.  Am.  J.  Phys.,  53:859. 

^  Yaghjian,  A.D.  (1980)  Electric  dyadic  Green's  functions  in  the  source  region.  Proc.  IEEE. 
68:258. 
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T - n- 

r-r' 


r  niHy^k  I  t-t'  I )  dx'dy'. 

Ja-a^ 


The  position  vector  r  has  been  replaced  by  the  two-dimensional  position  v^ector  t  =  ^  +  y y  . 
and  the  principal  volume  has  been  replaced  by  its  two-dimensional  analog,  the  cross-scctional 
principal  area  Aj.  Substituting  this  Hankel  function  expression  for  ij)  in  Eq.  (2).  we  have: 


E  (t) 


-1  /.  L.-J(t) 

- -  lim  f  [V.'xV.'x  J(t')]H„  (klt-t’l)  dA'  +  -  . 

4"»'„  .^-oJa-a, 


13) 


It  is  legitimate  to  move  the  double  curl  outside  the  z’-integration  so  it  only  includes  y.  because 
the  current  is  constant  with  respect  to  z'.  The  principal  volume,  of  course,  is  also  constant  in 
z  .  Thus,  the  limits  are  not  variable,  the  z'-integration  is  continuous,  and  the  derivative 
operator  can  be  moved  in  and  out  of  the  z'  integral  sign.  From  now  on  we  will  denote  1 1  -  t'  | 
as  p.  Using  a  vector  identity,  we  recast  curl  (curl  J)  and  write 


E  (t)  =  lim  f  J(t’)  •  (V,V,- vf)H‘„'\kp)  dA'  + 

.wo  J^-\ 


•  J(t) 

i">‘'o 


(4) 


With  some  algebra,  we  arrive  at  the  two-dimensional  scattered  field  expression 


-k  /» 

E  (t)  -  : -  lim  f  J(t')  ■  g"  dA’ 


I,  •  J(t) 

2i(i>FQ 


(5) 


where 


KZiZi 

4  I  k"^ 


+ 


H 


(1) 


(6) 
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and 


L 


6  " 


2  ■ 


(7) 


We  have  inserted  the  value  of  the  source  dyadic  lEq.  (7)1  corresponding  to  a  square  principal 
area  because  we  wish  to  divide  the  target  into  square  patches  later  in  the  numerical 

calculation.  If  a  different  patch  geometry  is  desired.  must  be  adjusted  accordingly.^ 

We  now  consider  the  two  possible  polarizations  of  the  incident  wave.  We  begin  with  the 
simpler  TM  polarization,  in  which  E  =  zE^.  For  a  plane  wave,  this  corresponds  to  either 
E  =  ze**^  or  E  =  ze'*^^.  We  need  only  solve  the  z-component  of  the  equation: 


-K  y 

E.(t)  =  -  lim  r  (J{t')  •  g"),  dA'. 

6^Ja-a, 


(8) 


The  equation  for  this  polarization  is  simplified  by  the  vanishing  of  the  source  dyadic  term.  In 
addition,  since  the  two-dimensional  derivative  operator  has  no  z-comi)onent  and  the  z- 
component  of  1  is  1.  the  Green  s  function  is  also  simplified: 


(9) 


The  equation  solved  by  the  FORTRAN  program  is  then 


Es, If”  r  dA'. 

ft_0  Ja-a., 
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To  evaluate  the  electric  field  numerically,  we  approximate  the  integral  as  a  summation 
over  the  fa( f  of  the  target,  livided  into  N  square  patches,  each  with  a  so\uee  i)oint  at  its  c  t  nter 
To  exclude  the  i)i  inci[)al  area,  we  initially  allow  it  to  occupy  an  entire  s(piare  pate  li.  and 
.simj)ly  omit  the  contribution  from  that  patch.  The  limit  occurs  as  the  number  of  patches 
increases  for  a  target  of  fixed  size.  In  Section  5  we  improve  the  accuracy  oi  tliis  treatment  to 
let  the  principal  area  approach  zero  while  the  square  patch  size  remains  fixed,  and  we  then 
calculate  the  contribution  from  the  area  between  these  two  boundaries.  In  the  meantime. 


G 


ilowever.  we  exclude  the  singular  point  by  omitting  the  entire  self-cell  from  tlie  siunmalion  as 
follows: 


AA(i)|iQ 


(k  1 1„  -t,  I  ;  i  xn  . 


(11) 


Here.  AA  is  the  patch  area,  constant  for  all  patches,  and  Jz  is  considered  constant  over  the 
patch,  at  its  center  value.  Since  we  want  to  consider  dielectric  targets,  having  only 
polarization  current,  we  substitute 


d|i)l  -  dliiil - i"lP  (12) 

=  -i(.i(K  -  r„)E,„, 

= 


Assuming  t  is  constant  over  all  patches,  our  scattered  field  expression  then  reduces  to 


AA(riH„T 

4 


(13) 


and  a  summation  equation  is  written  as  follows  for  each  observation  patch,  n: 


^lIK  Z  ^^11  ^ 


AA(i)H„T 


g  H'„'\klt, 


)E,„,  ,(l|i  ;  i  -n. 


114) 


From  a  physical  viewpoint,  if  we  were  to  stand  at  the  center  point  of  patch  n  and  measure  the 
field  contributions  from  all  the  other  patches  except  n  we  would  sense  the  following 
contributions: 
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AA{o^„t  ,  , 

n=  1  :  Ei,^  ,  (t, )  =  E,„,  ,  (t,)  +  - ^  [rf”(k  1 1,  -ta  I  )Et„,  ^it^) 

+  -t,|)E„„  ,(t.,) 

...+  Hy’(k|t, -tNl)E,„,  .(V) 
AAa)u„T  ,  , 

n=2:  E,,^.^(t2)=  E,„,  z  (ta)  +  - 5 -  tH<”(k|ta -t,  l)E,„,  ,(t,) 

+  H<„'\k|t2-t3l)E,„j,(t3) 
...+  H'“(k|ta  -tNl)E,o,  ,{tN)]. 


Clearly,  we  may  substitute  unity  for  the  self-cell  position  (that  is.  the  diagonal  matrix  element)  in 
our  system  of  equations  and  we  can  write  it®  in  the  form  Ax  =  B.  For  convenience,  we  then 
multiply  Eq.  ( 1 5)  by  4/AAcop,,.  and  replace  E,„|  by  Jz/t.  The  N  x  N  matrix  equation  appears  in  the 
code  in  terms  of  J,,.  with  the  constant  4/AA{opot  on  the  matrix  diagonal: 


4.  THEORY  FOR  TE  POLARIZATION 


The  TE  incident  wave  is  polarized  in  the  perpendicular  direction,  so  E=  x  E^+  y  E^,  and  the 
source  dyadic  does  not  vanish.  In  addition,  the  integral  involves  coupling  between  the  x-and  y- 


®  Su,  C.C.  (1987)  Calculation  of  electromagnetic  scattering  from  a  dielectric  cylinder  using 
the  conjugate  gradient  method  and  FFT,  IEEE  Transactions  on  Antennas  and  Propagation, 
AP-35:1418. 
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directions,  whereas  tlie  TM  polarization  was  completely  decoupled.  Our  scattered  field 
equation  is 


i(r)P(j  liin  f  T  f  *  ~V")  ■  I  *■“*'  1 

fe-ojA-Aj**  V  k  / 


J,(t) 

2i(p)t- 


(17) 


where  the  dot  product  in  the  integrand,  upon  manipulation,  becomes 


fi  + 


V  V  \ 
I  t 


J 


X 


1  I  i)  uO 

<'x  !)y  ° 


-L  J 
k^  " 


A  Ah^‘^ 

ay  ax  0 


+  j 

y 


(18) 


We  refer  to  the  RHS  of  Eq.  ( 18)  as  Vt-  Substituting  the  scattered  field  expression  into  the 
standard  equation,  we  have 


'lot 


IN 

i(‘)  =  >‘m  )  V,  [p.  E,  AA„ 

AA-*0 


2iii)f  n 


(19) 


The  last  term  is  the  source  dyadic  corresponding  to  the  square  principal  area,  as  before. 
Rearranging,  we  have 


^Inc  t^*n^  ”  ^tol  1^*11' 


c(t. 


2i(i)ff 


wPqAA 


N 

I"''" 


't  lot 


(t')]  ;  i;<n 


(20) 


for  each  of  N  patches.  Again,  the  stipulation  i  x  n  excludes  the  self-patch,  where  n  denotes  the 
patch  containing  the  observation  point  and  i.  the  patch  with  the  source  point.  The 
components  of  Eq.  (20)  are: 
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Inc^^n^  ~ 


1 

T(t„) 


Ey  =  Jy(tn) 


c(t.. 


1 

2i(U£Q 


+ 


wMqAA 

4 


+  J 


1  _2 _ Lh"^ 

y'i^  ox  ay  “ 


(21) 


1 

2i(o£Q 


idUoAA 

4 


+  J 


1  a  a 
ax  ay 


Note  that  the  argument  of  the  Hankel  function  is  always  kp  Indicating  that  each  component  of 
the  individual  electric  field  Induces  both  components  of  current.  J  —  the  reason  why  the  x-and 

y-components  remain  coupled.  Difierentlating  the  Hankel  function.  h|".  we  can  express  the 
summation  in  the  Ex  mr  equation  of  Eqs.  (21).  as: 


H, 


H2  (  x„  -  X,) 


+J., 


H2 


(Xn-  x,)(yn  -  y,) 


(22) 


A  similar  result  occurs  for  the  summation  in  the  Ey  m,.  equation.  If  we  then  write  out  the 
summation  for  a  target  whose  cross-section  is  divided  into  N  square  patches,  we  obtain  this 
time,  a  2N  x  2N  matrix  equation,  whereas  the  TM  polarization  had  only  an  N  x  N  matrix  to 
invert. 
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We  do  not  consider  oblique  incidence,  but  only  the  two  cases  =  ye*''’^  (vertical  target)  and 
=  xe**'^  (horizontal  target).  The  constant  C  on  the  matrix  diagonal  is 


fi  -  1  4 

[x  2iwef,  AAcoiij^ 


(24) 


After  inverting  the  matrix,  we  then  add  Jxh  and  to  get  (p’).  To  calculate  the  far-field  RCS,  we 
use  the  two-dimensional  formula. 


.  r**. , , 

v '  11  la  > 


(25) 


where  the  incident  H  field  in  free  space  is  Hjik-  =  -(Eo/Mo)*'^^  Einc  and 


)  AA 


(26) 


II 


Thus,  our  final  equation  for  the  scattering  cross  section  is 


^2-0  ^*1’^ ' 


/  AA\2  1  .  V 

4^(^) 


J,e 


lkcos(^  -  ♦,)  p, 


(27) 


5.  DERIVATION  OF  SELF-CELL  CORRECTION 


Until  now.  we  have  excluded  the  entire  self-cell  as  though  it  were  the  limiting  principal 
area  when,  in  fact,  it  may  be  too  large  to  accurately  model  a  tiny  area  shrinking  in  toward  the 
source  point  at  its  center.  The  accuracy  improves  as  the  number  of  cells  per  wavelength 
increases,  or  the  relative  cell  size  shrinks.  However,  this  is  more  and  more  computationally 
demanding.  As  an  alternative  method  for  improving  the  accuracy  of  our  calculation,  let  us 
consider  first  the  analogous  three-dimensional  situation,  involving  the  volumetric  self-cell. 
Vgc.  surrounding  a  principal  volume.  Vj,  (see  Figure  2).  Since  we  wish  to  use  V;,  instead  of  V^, 
for  the  excluded  region,  we  must  calculate  the  contribution  to  the  total  field  from  the  region 
between  Vgc  and  Vj,. 
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For  these  calculations,  we  will  use  spheres  and  circles  rather  than  cubes  and  squares; 
conveniently,  both  shapes  have  identical  principal  volume  and  principal  area  expressions.' 
The  equation  we  have  used  until  now  for  the  scattered  E  field  excluding  an  entire  self  cell  in 
place  of  V()  is 


XLj,  •  J(r) 
G[.-  (r.r  )  •  J(r')  dV  +  — - 


(28) 


We  prefer  to  use 


E(r) 


iii)|.i,,  r  G  •  J  dV 

J\-\- 


(29) 


where  the  first  term  represents  the  contribution  from  all  the  cells,  the  second  term  repre.sents 
the  principal  volume,  and  the  third  term  represents  the  'intermediate"  contribution  we  will 
proceed  to  calculate. 

Following  the  ideas  outlined  by  Nachampkin.^  we  formulate  AE|r)  as: 


AE(r)  =  lim  imu,,  C  G,.  (r,  r')  •  J(r')  dV. 


(30) 


Assuming  J(r')  is  constant  over  Vg,.. 


AE(r)  =  lim  imu,,  J(r')  T  G,,  (r.  r')  dV 


(31) 


where 


^  Nachampkin.  J.  (1989)  Integrating  the  dyadic  Green's  function  near  sources.  IEEE 
Transactions  on  Antennas  and  Propagation.  38:919. 
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G 


(32) 


Vx  Vx(\i»I) 

4jii? 


We  may  easily  express  the  electric  Green's  function  in  terms  of  the  magnetic  Green's  function: 


G 


nj 


vxt^i 

4x 


(V  X  i;.!). 
4x 


(33) 


Using  Stokes'  theorem,  we  then  convert  the  volume  integrals  into  surface  integrals: 


AE(r)  =  lim  — ^ —  J 
IwEn  ■ 


Eliminating  G^  from  the  integrand  via 


n  X 


G„,(r.  rO  dS*  - 


X." 


G„,(r.  r  )  dS' 


(34) 


iV  X  G, 


4x 


(iV  X  X  i) 


(35) 


with 


(36) 


gives 


n'  X  G„,  = 


(ikU  -  lie 
(4nR)^ 


lkl< 


[n'R  -(n'  R)  I]  . 


(37) 


Next,  we  integrate  this  quantity  over  the  surfaces  S^,  and  S*.  If  we  choose  V’^,,  and  V*  to  lie 
spheres,  where  the  radius  of  V^,  =  a,  and  that  ofVj  =  rj.  when  we  integrate  over  the  surface.  S^,  . 
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the  quantity  R  =  Ir  -  r  =  a.  since  the  observation  point  Is  at  the  center  of  the  sphere.  V^, .  and 
the  source  points,  r'  are  on  Its  surface  (see  Figure  3).  Likewise,  when  we  Integrate  over  the 
surface  Sg.  R  =  rg. 


ft 


^SC 


Figure  3.  Principal  Volume  Relative  to  Volume  of  Self-Cell 


I 


n'  X  G^(r.  r  )  dS'  =  (tka-1)  e'**^ 

4n 


dS' 


-f  ^ 

471  Js  a^ 

n' 


(38) 


We  recognize  the  term 


(391 


as  the  definition  of  and  since  I  Is  constant,  the  second  Integral  reduces  to  I  ,  Thus. 


1 


h  X  G^tr.  r  )  dS'  =  —  (lka-1)  e'**"  [-L„  +  I] 

m  sc 


(40) 
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Analogously,  with  rg  as  the  radius  of  Vg. 


1 


n'  X  G„(r.  r  )  dS'  =  ^  (IkTs-l)  e"^«  l-L^  +  II. 
47t 


(41)  • 


So.  our  equation  for  the  Intermediate  electric  field  becomes 


AE(r)  =  llm 
6-»o 


4jtla)eo 


{  (lka-1)  6“^“  [-L,^  +  I] 


-  (IkTg-l)  e“^»  l-Lj  +  I]  }. 


(42) 


Taking  the  limit  as  8-)0.  so  rg-)0.  and  considering  the  source  dyadic  for  a  sphere. 


L 


8 


I 

3  ■ 


(43) 


We  now  have: 


AEir)  =  — {  (Ika  +  1}. 

6rcicae„ 


(44) 


Now.  we  solve  the  two-dimensional  problem  In  a  completely  analogous  manner,  using  a 
circular  self-cell  and  circular  principal  area.  For  the  TM  polarization,  we  had  the  equation  for 
the  scattered  electric  field. 


E^(t)  =  lean 


J 


(g“  (t.  t')  •  J  (t')U  dA’  + 


AE^tt). 


(45) 


To  develop  an  expression  for  the  intermediate  field,  we  once  again  assume  the  current  Is 
constant  over  the  area  of  the  self-cell. 
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(■10) 


AE^d)  =  lirn  ioiHg  f  dA'. 

6-0  JavA 


and  express  the  Green's  function  in  terms  of  tlie  Hankel  function  in  the  integrand: 


-(DU  ^  . 

AE^(t)  =  — ^  lim  r  H„Up)  dA'. 

^  5^0 


(•17) 


We  can  solve  the  indefinite  Integral: 


(48) 


If  we  put  in  the  limits  and  a,  we  have  our  result  -- 


AE,  (t)  =  J,  llm[aHV’(l<a)  -  pjH'J^kpj)!. 

6—0 


(49) 


Continuing  with  the  TE  polarization,  we  had  the  equation  for  the  scattered  field: 


E(t)  =  i<.)p„  f  g°  (t.  t’)  •  J  (f)  dA' 
Ja-a^ 


J,(t) 

2i<.)E„ 


+  AE{t). 


(50) 


and 


AE{t)  =  idipo  lim  f 

d-*o  Ja, 


I 


(V, 


X  V,  X 


•J,  dA'. 


(51) 


just  as  in  the  three-dimensional  case.  Next,  we  write  AE  in  terms  of  g",.  assuming  J,  is 
constant  over  A^, 
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AE(t)=— (Vtxg°)  dA'  -  lim  f  (VtXg°)  dA 
Ja_,  8-»0 


Since  we  can  also  show 


f  (Vj  X  A)  dA  =  [  (fi,  X  A)  dC. 
Ja  Jc 


we  apply  this  two-dimensional  version  of  Stokes'  theorem  to  the  two  area  Integrals  In  Eq.  (52), 
to  get 


=  *  [  f  (n'  ><  g  m)  dC  -  llm  f  (h'  X  g  ° 


)  dC  . 


Substituting  Hankel  functions  for  the  cross  product  In  this  case  gives 


n'  X  g°  =  4  h'o'  -  il. 


The  gradient  of  the  Hankel  function  comes  out  to  be 


T7L.<i»  _  (t  -t’) 

TTTT  ‘  ■ 


n,  X  g  =  -L-  H*/’  li 


and  substituting  back  Into  the  contour  Integrals, 
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AE(t) 


4(iif 


h'i” 


X. 


(-n,e,  +  i)  dC 


-  limH'/’r 
6-0 


(-n,e,  +  I)  dC 


(58) 


Note  tliat  the  Hankel  function  is  constant  over  tiie  perimeter  contours,  as  we  evaluate  tlie 
four  Integrals  around  the  circles  Csc  and  Cj,.  with  the  radius  of  Csi-  =  a,  and  that  of  C,s  =  (i.s.  I  is 
also  constant,  and  we  recognize 


n,e, 

2.iT' 


dC 


h- 


(5f(l 


For  a  circle,  we  use^ 


h  =  Isr 


I 

2' 


(GO) 


and  our  result  becomes 


AE(t)=-^J. 


aH‘i’(ka)-  lim  gj  h',"  (kp^) 
6-0 


(61) 


This  is  quite  similar  to  the  result  for  the  TM  polarization.  For  the  last  term,  we  can  use  the 
small  argument  approximation  for  the  Neumann  function  to  get 


lim  -  ip,^  1  r(l )  k(),.,) 

6-0  ^ 


(62) 


The  gamma  function  evaluates  to  unity,  giving  us  the  value 


2i 

-Tk  (63) 
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for  the  desired  limit.  Substituting  this  into  our  results  for  both  two-dimensional  polarizations 
we  have; 


AE,(t) 


-(OPoJl 

2k 


aH*!*  (ka)  + 

ak 


(64) 


and 


AE(t)  = 


-(OPoJl 

4k 


aH'/'lka)  +  -21 
‘  jtk 


(65) 


In  the  numerical  calculation  we  simply  replace  the  circular  principal  area  with  an  equivalent 
square  principal  area  by  forcing  s  =  or  in  the  case  of  volumes  in  three  dimensions, 

s  =  (4/3  .'i)*''^a. 


6.  NUMERICAL  CALCULATION  OF  SCATTERING  CROSS  SECTION 


Once  we  have  arrived  at  the  matrix  equation,  we  simply  choose  a  matrix  solution  method 
to  calculate  the  vector  x  in  Ax  =  B,  corresponding  to  the  current  induced  in  the  target.  Recall 
also  that  B  is  the  incident  field  vector,  and  the  matrix  A  consists  of  Hankel  functions  that 
depend  only  on  the  distances  between  the  source  and  observation  points.  We  initially  chose  a 
Gauss-Seidel  routine  to  invert  the  matrix,  but  recommend  the  conjugate  gradient  method  be 
used  in  a  continuation  of  this  project.®  The  program  discussed  in  this  report  works  only  for  a 
rectangular  cylinder  divided  evenly  into  square  patches,  but  can  be  modified  for  alternative 
geometries;  the  analytical  equations  hold  for  an  arbitrarily-shaped  scatterer.  It  might  be  a 
useful  exercise  to  modify  the  code  so  it  will  accommodate  a  circular  cylinder,  since  a  closed 
form  solution  exists  for  comparison  with  that  particular  case,  although  the  physical  optics 
standard  is  thought  to  be  a  sufficiently  accurate  basis  for  comparison. 

Figure  4  shows  the  model  on  which  the  code  is  based.  A  plane  wave  (at  left)  is  incident  on 
the  infinite  rectangular  cylinder,  and  the  bistatic  RCS  is  calculated  as  a  function  of  the  angle  ^ 
as  indicated.  The  cross-sectional  area  of  the  cylinder  is  divided  evenly  into  square  patches, 
indexed  according  to  Figure  5.  to  take  advantage  of  the  rectangular  symmetry.  The  origin  of 
the  coordinate  system  is  always  set  up  to  be  at  the  geometric  center  of  the  target.  The  target 
dimensions  are  specified  on  input  and  the  program  adjusts  the  incident  field  vector 
accordingly  if  the  polarization  is  specified  as  TE.  The  program  expects  eitlier  a  square  or  a 


20 


2 


Figure  4.  RCS  Model  Used  in  Code 


vertical  target  that  Ls,  with  the  longest  side  in  the  y-direction.  The  z-diinension  is  infinite  for 
a  two-dimensional  problem.  The  incident  wave  vector  points  in  the  +x-direction.  The 
program  checks  to  make  sure  the  user  specified  a  vertical  target.  If  the  user  inputs  a  hori¬ 
zontal  target,  the  program  rotates  the  target  and  sets  the  incident  wave  vector  in  the  +3'- 
direction,  for  grazing  incidence.  Thus,  the  user  defines  the  two  cases  of  plane-wave  incidence 
by  the  order  in  which  the  dimensions  of  the  target  appear  in  the  input  file. 


Figure  5.  Indexing  Scheme  for  a  Rectangular  Target  Divided  Evenly  Into  Square  Patches 


Results  were  obtained  for  a  variety  of  targets.  Initial  test  runs  were  made  with  a  square 
cylinder.  0.8  x  0.8  wavelength,  for  z-polarization.  for  several  values  of  the  ratio  r/io.  When 
this  ratio  was  set  to  3.0.  runs  were  made  for  N  =  4.  16.  64.  and  100.  Results  for  these  runs 
aiipear  in  Figures  6  and  7.  The  ratio  was  then  set  to  1.1  to  simulate  a  target  that  was  almost 
invisible  to  the  wave  (free-space  target),  to  make  sure  the  RCS  approached  zero  Figure  8 
verifies  this  limit.  The  run  time  for  a  target  with  100  patches  (N  =  100)  was  tibout  10  minutes. 
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RCS  OF  AN  INFINITE  SQUARE 
CYLINDER  .8  x  .8  LAMBDA.  e/Eo  =  3,  ZPOL 
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Figure  6.  RCS  of  an  Infinite  Square  Cylinder.  0.8  x  0.8  k,  t/to  = 
Z  Polarization,  N  =  4.  16,  and  64 


RCS  OF  AN  INFINITE  SQUARE 
CYLINDER  .8  x  .8  LAMBDA,  e/Eq  =  3,  ZPOL 


(SOd  IVlOl)  Vl^OIS 
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Figure  7.  RCS  of  an  Infinite  Square  Cylinder.  0.8  x  0.8  X.  e/eq  = 
Z  Polarization.  N  =  4.  64.  and  100 


RCS  OF  AN  INFINITE  SQUARE 
CYLINDER  .8  x  .8  LAMBDA.  e/£o  =  1  -1 .  ZPOL 


Figure  8.  RCS  of  an  Infinite  Square  Cylinder.  0.8  x  0.8  k.  e/tQ  =  1.1 
Z  Polarization 


The  next  set  of  runs,  at  the  top  of  Figure  9,  also  had  a  z-polarlzed  Incident  wave,  but  the 
parameters  were  set  to  correspond  to  Richmond's  example  (bottom).  Richmond^  used  an 
epsilon  ratio  of  4.0  and  a  slab  measuring  50  patches  x  1  and  2.5  x  0.05  wavelengths.  We  used  a 
similar  slab  measuring  100  patches  x  2,  but  the  real  dimensions  are  still  2.5  x  0.05 
wavelengths.  Since  we  have  4  times  as  many  patches,  we  expect  an  Improvement  In  accuracy'. 
We  ran  both  orientations  with  respect  to  the  Incident  wave,  as  he  did.  We  refer  to  these  as 
"horizontal  target  "  (grazing  Incidence)  and  "vertical  target",  (normal  Incidence)  and  they  are 
depicted  In  the  two  columns  of  Figure  9.  (Although  we  use  (J»  =  0  as  backscatter  and  180"  as  the 
forward  scatter  direction  our  angles  In  Figures  9  and  10  have  been  shifted  to  conform  to 
Richmond's  definition.)  Figure  10  shows  the  corresponding  results  for  the  xy-  or  TE  polari¬ 
zation.  All  the  results  of  Figures  9  and  10  Include  the  self-cell  correction. 

The  graphs  are  Identical  for  z-polarlzatlon,  with  a  horizontal  target.  For  both 
polarizations,  the  vertical  target  graph  Is  closer  to  the  physical  optics  approximation  than  Is 
Richmond's.  The  horizontal  target  RCS  for  the  jty-polarizatlon  is  slightly  higher  than 
Richmond's,  but  he  gives  no  physical  optics  result  for  comparison  In  this  case. 

We  experimented  with  leaving  out  the  L*  J  source  term  for  the  TE  case,  and  found  that  as 
expected,  the  L  •  J  term  Is  necessary  to  achieve  a  correct  result.  The  self  cell  correction  then 
fine-tunes  the  accuracy  of  the  RCS.  We  have  presented  a  complete  development  of  the  theory' 
leading  to  this  calculation,  where  Richmond  omits  this  Issue  from  his  discussion,  or 
Incorporates  the  Information  contained  In  these  terms  by  some  other  means.  All  the  data 
were  produced  using  either  the  VAX  1 1-750  at  RADC/EECT  or  the  VAX  8650  scientific  machine 
at  AFGL. 


7.  FUTURE  EXTENSIONS 


The  original  Intent  of  this  project  was  to  continue,  and  to  solve  the  problem  in  three 
dimensions  as  well,  and  to  allow  for  an  Inhomogeneous  target  (where  n  and  E  vary  with 
position).  The  three-dimensional  problem  will  involve  Inverting  a  matrix  with  N  proportional 
lo  (d/X)'^  rather  than  (d/X)^  for  two  dimensions.  Although  It  has  been  said  that  the  anaMlcal 
calculation  Is  less  cumbersome  than  what  we  have  already  done  for  the  two-dimensional  case, 
the  computer  resources  needed  are  significantly  greater. 

Following  are  some  suggestions  for  decreasing  the  computing  time.  First,  one  should  use  a 
more  efficient  matrix  solver,  (such  as  conjugate  gradient  Iteration  combined  with  the  FFT). 
and  consider  the  advantage  of  using  dimensionless  variables  In  the  code,  which  also  simplifies 
the  process  of  varying  parameters.  Furthermore,  this  problem,  especially  the  Inhomogenous 
form.  Is  a  good  candidate  for  a  parallel  processor.  It  might  also  be  Interesting  lo  consider 
doing  a  tlmestepped  simulation  of  wave  scattering,  as  an  alternative  to  using  the  computer  In 
Its  usual  number-crunching  role. 


26 


Figure  9.  Upper  Lelt:  RCS  of  an  Infinite  Square  Cylinder.  0.05  x  2.50  r/fQ  =  4.0.  Z  Polari¬ 
zation.  With  Vertical  Target.  Upper  right:  with  horizontal  target;  for  comparison  with  upper 
left.  Lower  left:  calculated  scattering  patterns  of  a  homogeneous  plane  dielectric  slab  of  the 
same  thickness  and  width  as  above,  with  a  plane  wave  having  normal  incidence,  and  lower 
right:  with  a  plane  wave  at  grazing  incidence  (after  Rlchmondl 
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